Catalyst-like modulation of transition states for CFTR channel opening 
and closing: New stimulation strategy exploits nonequilibrium gating 
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Cystic fibrosis transmembrane conductance regulator (CFTR) is the chloride ion channel mutated in cystic fibrosis 
(CF) patients. It is an ATP-binding cassette protein, and its resulting cyclic nonequilibrium gating mechanism sets 
it apart from most other ion channels. The most common CF mutation (AF508) impairs folding of CFTR but also 
channel gating, reducing open probability (P„) . This gating defect must be addressed to effectively treat CF. Com- 
bining single-channel and macroscopic current measurements in inside-out patches, we show here that the two 
effects of 5-nitro-2-(3-phenylpropylamino)benzoate (NPPB) on CFTR, pore block and gating stimulation, are 
independent, suggesting action at distinct sites. Furthermore, detailed kinetic analysis revealed that NPPB potently 
increases P t „ also of AF508 CFTR, by affecting the stability of gating transition states. This finding is unexpected, 
because for most ion channels, which gate at equilibrium, altering transition-state stabilities has no effect on P Q ; 
rather, agonists usually stimulate by stabilizing open states. Our results highlight how for CFTR, because of its 
unique cyclic mechanism, gating transition states determine P„ and offer strategic targets for potentiator com- 
pounds to achieve maximal efficacy. 



INTRODUCTION 

CFTR is a chloride ion channel (Bear et al., 1992) crucial 
for the salt water balance of several polarized epithelia. 
Mutations in CFTR are the cause of cystic fibrosis (CF) 
(Riordan et al., 1989), the most common lethal genetic 
disease among Caucasians, an incurable, devastating 
multi-organ disorder (O 'Sullivan and Freedman, 2009). 
The most common CF-causing mutation (carried by >90% 
of patients), deletion of phenylalanine 508 (AF508), 
severely impairs both protein folding/stability (Cheng 
et al., 1990; Okiyoneda et al., 2010; Rabeh et al., 2012) 
and channel open probability (P Q ) (Wang et al., 2005; 
Miki et al., 2010). Thus, much effort is focused on the 
development of "correctors," chemical chaperones that 
promote folding/stability of the AF508 CFTR protein, 
and of "potentiators," which stimulate P Q of AF508 (or 
other mutant) CFTR channels. 

CFTR is an ATP-binding cassette (ABC) protein, with 
a typical ABC architecture consisting of two trans- 
membrane domains (TMDs) and two cytosolic nucleo- 
tide-binding domains (NBDs) . The additional cytosolic 
regulatory (R) domain, unique to CFTR (Riordan et al., 
1989), is a substrate for cyclic AMP-dependent PKA, 
and its phosphorylation is required for CFTR channel 
activity (Tabcharani et al., 1991). CFTR's chloride ion 
pore is formed by the TMDs and gated by a cycle of 
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ATP binding and hydrolysis at the two NBDs (Gadsby 
etal., 2006). 

Despite the functional divergence of the ion channel 
CFTR from other ABC proteins, which are mainly active 
transporters, a conserved molecular mechanism cou- 
ples cycles of conformational changes in NBDs and in 
TMDs of all ABC proteins (Locher, 2009). In the pres- 
ence of ATP, NBDs of ABC proteins form tight head- 
to-tail dimers, occluding two molecules of ATP; each 
nucleotide is sandwiched between the conserved Walker 
A and B motifs of one NBD and the conserved signature 
sequence of the other, which together form a compos- 
ite catalytic site with ATPase activity (e.g., Karpowich 
et al., 2001; Smith et al., 2002). NBD dimers glued to- 
gether by two ATP molecules are extremely stable, but 
they dissociate after ATP hydrolysis (Moody et al., 2002; 
Verdon et al., 2003) . In ABC-C family members, includ- 
ing CFTR (ABCC7) , the composite binding site formed 
by Walker motifs of NBD1 and signature sequence of 
NBD2 ("site 1") is degenerate and catalytically inactive 
(Aleksandrov et al., 2002; Basso et al., 2003); ATP hy- 
drolysis occurs only in "site 2" (formed by Walker motifs 
of NBD2 plus signature sequence of NBD1 ) . In full-length 
ABC exporter structures, tightly associated NBD dimers 
are linked with outward-facing (Dawson and Locher, 
2006; Ward et al., 2007) and fully or partially separated 
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NBDs with inward-facing (Ward et al., 2007; Aller et al., 
2009; Hohl et al, 2012) TMD conformations, suggest- 
ing that formation of ATP-bound NBD dimers flips the 
TMDs from inward to outward facing, whereas ATP hy- 
drolysis drives dimer dissociation to reset the TMDs to 
inward facing (Locher, 2009) . 

In the ion channel CFTR, the NBD dimerization-dis- 
sociation cycle is coupled to the opening and closure of 
the transmembrane permeation pathway: the pore opens 
to a burst when the tight NBD dimer forms, and it closes 
from a burst when the dimer interface separates around 
site 2, after hydrolysis of the ATP bound there (Vergani 
et al., 2005). (CFTR current records display a bursting 
pattern: open events, which last for ^2-300 ms, are inter- 
rupted by brief [<10-ms] ATP-independent "flickery" clo- 
sures, distinct from the long closed events [^1 s] that 
reflect NBD dimer separation ["interburst closures"; 
Vergani et al., 2003]; in this study, the terms "pore open- 
ing" and "pore closure" are used synonymously with en- 
tering and exiting a burst.) In the absence of hydrolysis, 
pore closure is extremely slow: site-2 mutations that dis- 
rupt ATP hydrolysis slow the closing rate by ^100-fold 
(Gunderson and Kopito, 1995). Thus, gating of WT 
channels is an essentially unidirectional cycle; channels 
that open to a prehydrolytic open state (Oi) preferen- 
tially progress to a post-hydrolytic open state (Oo), to 
then close through a pathway (G\— >Cv) distinct from 
pore opening (Q— >Oi) (shown in insets throughout 
Figs. 3-5; Csanady et al., 2010) . This far-from-equilibrium 
operation is rare among ion channels but is an essential 
property of ABC transporters that mediate unidirectional 
uphill transport. At saturating [ATP], the CFTR func- 
tional cycle contains two rate-limiting steps with relatively 
high free energy barriers: step Ci— >Oi (rate of ^T s _1 ) 
determines opening rate, whereas step Oi— >0% (rate of 
^4 s _1 ) rate limits closure (Csanady et al., 2010). 

Because the TMDs are the most divergent parts of ABC 
proteins, they are promising targets for the development 
of highly selective drugs. The arylaminobenzoate 5-nitro- 
2-(3-phenylpropylamino)benzoate (NPPB) is an anionic 
compound that binds to CFTR's TMDs, and it has been 
studied extensively for its pore blocker properties (Zhang 
et al, 2000; Zhou et al., 2010) . Unexpectedly, Wang et al. 
(2005) discovered that NPPB, in addition to blocking the 
pore, dramatically increases P Q (i.e., potentiates gating) of 
both WT and AF508 CFTR Given the potential therapeu- 
tic impact of CFTR potentiators, the present study aimed 
to identify the exact molecular mechanism by which 
NPPB stimulates P„. As a control, we compared the effects 
of NPPB to those of MOPS, another well-characterized 
CFTR pore blocker (Ishihara and Welsh, 1997), but with- 
out potentiating effects. Our results offer a conceptually 
novel approach to robustly stimulate CFTR that exploits 
its unique nonequilibrium gating cycle, and that differs 
fundamentally from strategies applied so far to modulate 
the activity of any ion channel. 



MATERIALS AND METHODS 

Molecular biology 

WT and mutant CFTR cDNA in pGEMHE was transcribed in vitro 
using T7 polymerase, and 0.1-10 ng cRNA was injected into Xeno- 
pus laeuis oocytes as described previously (Csanady et al., 2010). 

Electrophysiology 

Current recordings were done at 25°C in inside-out patches ex- 
cised from oocytes 2-3 d after RNA injection. Pipette solution 
contained (mM): 136 NMDG-C1, 2 MgCl 2 , and 5 HEPES, pH 7.4 
with NMDG. Bath solution (mM: 134 NMDG-C1, 2 MgCl 2 , 
5 HEPES, and 0.5 EGTA, pH 7.1 with NMDG) was continuously 
flowing and could be exchanged with a time constant of 20-30 
ms. CFTR channels were prephosphorylated by 1-2-min exposure 
to 300 nM PKA catalytic subunit (Sigma-Aldrich). Recordings 
were done in the presence of saturating (2 mM) MgATP; for the 
K1250A mutant 10 mM MgATP was used to compensate for its 
greatly impaired apparent ATP affinity (Vergani et al., 2003). 
MgATP (Sigma-Aldrich) was added from a 400-mM aqueous stock 
solution (pH 7.1 with NMDG), and Na 4 -pyrophosphate was added 
from a 500-mM aqueous stock solution, together with equimolar 
MgCl 2 . NPPB (Sigma-Aldrich) was added from a 250-mM stock 
solution in DMSO; by spectrophotometry, NPPB solubility in 
water was >400 pM. Solutions containing MOPS were adjusted to 
pH 7.2; thus, [MOPS~] was ~50% of total [MOPS] (pK„ = 7.2). 
Currents recorded (Axopatch 200B; Molecular Devices) at a 
bandwidth of 2 kHz were digitized at 10 kHz (Digidata 1322A; 
Pclamp9; Molecular Devices) . 

Data analysis 

For gating analysis, current records were refiltered at 100 Hz. 
Steady-state mean burst and interburst durations in patches con- 
taining fewer or equal to two channels (Fig. 6 F) were obtained by 
maximum likelihood fits of the closed-open-blocked scheme to 
the ensembles of all dwell-time histograms (Csanady, 2000) ; this 
approach efficiently separates brief ATP-independent "flickery" 
closures from long "interburst" closed events (Csanady et al., 
2010). For burst analysis of one-channel records (Fig. 5), brief 
closures were suppressed using the method ofjackson et al. ( 1983) . 
Distributions of burst durations longer than t lm . = 12 ms were 
fitted to gating models by maximum likelihood (Csanady et al., 
2010) ; alternative models were compared using the log-likelihood 
ratio test (Csanady, 2006). Average unitary current amplitudes 
were estimated from Gaussian fits to amplitude histograms of 
heavily (10 Hz) filtered records (see Fig. 1 B). 

Fitting of pH titrations 

Because titration of NPPB (Fig. 9 A) had to be done in a dilute 
(100-uM) solution, the titration curve was fitted by accounting for 
the presence of ~10 pM of aqueous C0 2 . Thus, the data were fit- 
ted to a curve expected for a mixture of two weak acids: for one 
component (C0 2 ) a fixed concentration of 10 pM and pK„* = 
6.47 was used, whereas concentration and pK„ of the second com- 
ponent (NPPB) was left free during the fit. The slight deviation of 
the data from the fitted curve at pH >9 (Fig. 9 A) reflects inevita- 
ble progressive accumulation of dissolved C0 2 at very basic pH. 

Statistics 

All symbols and bars represent the mean ± SEM of at least five 
measurements, obtained from several patches. For parameters 
extracted from macroscopic currents, the average number of 
measurements was greater than nine, and the average number of 
patches was greater than five, for each data point (bar or symbol) 
shown. For the single-channel data shown in Figs. 5 and 6 F, the 
numbers of patches analyzed for each condition are provided in 
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the figure legends. Statistical significance was evaluated using Stu- 
dent's two-tailed «test (*, P < 0.05; **, P < 0.01). 

Online supplemental material 

Fig. SI illustrates the kinetics of pore block by NPPB and MOPS. 
Fig. S2 demonstrates the effect of NPPB on the unlocking rate of 
WT CFTR channels locked open by ATP plus pyrophosphate. The 
online supplemental material is available at http://www.jgp.org/ 
cgi/content/full/jgp.201311089/DCl. 



RESULTS 

Voltage-dependent pore block by NPPB and MOPS 
reduces average unitary current through bursting 
CFTR channels 

Cytosolic anionic pore blockers are driven into the CFTR 
pore by hyperpolarizing voltages and disrupt the flow 
of chloride ions, causing flickery block. To deconvolve 
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Figure 1. Voltage and dose dependence of pore block by NPPB and MOPS - . (A and B) Unitary WT CFTR currents in symmetrical 
^140-mM CP with and without cytosolically applied blockers, displayed at two bandwidths (f„ filter corner frequency). We could not 
reliably discern two distinct conductance levels in MOPS" (Gunderson and Kopito, 1995). (C) Apparent unitary amplitudes at 10 Hz, 
with and without blockers, plotted against voltage. (D) Fractional unitary current in 210 uM NPPB and 20 mM MOPS - (symbols) and 
Boltzmann fits (solid lines) ; midpoint voltages (V1/2) and apparent valences (z) were V1/2 = — 85 ± 3 mV and z = —0.45 ± 0.03 for MOPS - 
and V1/2 = 24 ± 2 mV and z = -0.51 ± 0.03 for NPPB. (E, G, I, and K) Slowly decaying macroscopic "locked-open" currents of E1371S 
CFTR in the absence of bath ATP; brief applications of 210 pM NPPB or 20 mM MOPS - at various voltages (E and G) or of various 
blocker concentrations at —120 mV (I and K). Dotted lines show zero-current levels estimated from final segments; in E and G, small 
(<1 pA/40 mV) linear seal currents were subtracted. (F and H) Voltage dependence of macroscopic current block by (F) 210 pM NPPB 
and (H) 20 mM MOPS - (closed symbols) and Boltzmann fits (solid lines); V 1/2 = -87 ± 3 mV and z = -0.54 ± 0.07 for MOPS - and V 1/2 = 
7 ± 3 mV and z= —0.48 ± 0.03 for NPPB; open symbols and dotted lines replotted from panel D. (J and L) Dose dependence of macroscopic 
current block by (J) NPPB and (L) MOPS - at —120 mV (closed symbols) and Michaelis-Menten fits (solid lines). Open symbol in L: 
block by 40 mM of total MOPS (MOPS-H plus MOPS - ) at pH 6.2; calculated [MOPS ] = 3.6 mM: fractional block by MOPS depends on 
[MOPS - ], not total [MOPS], confirming anionic MOPS - to be responsible for pore block (Ishihara and Welsh, 1997). 
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effects on gating from those on permeation (pore block) , 
we first quantitated the latter in inside-out patches of 
single WT CFTR channels gating in 2 mM of bath ATP 
(Fig. 1 A, left). Bath application of 210 pM NPPB elicited 
flickery closures whose frequency increased at more neg- 
ative potentials to such an extent that distinct unitary gat- 
ing events disappeared, leaving "fuzzy" noise (Fig. 1 A, 
middle) . In contrast, bath application of 20 mM MOPS - 
appeared to reduce unitary current amplitudes at nega- 
tive voltages (Fig. 1 A, right). This seemingly different 
behavior reflects the lower affinity of MOPS - for the 
CFTR pore, resulting in briefer interactions and hence 
blocked events too fleeting to be resolved at our record- 
ing bandwidth of 100 Hz (Fig. SI; Ishihara and Welsh, 
1997; Zhang et al., 2000) . On the other hand, after filter- 
ing the data at 10 Hz (Fig. 1 B), the effect of NPPB on 
unitary currents appears similar to that of MOPS ; both 
appear to reduce unitary amplitudes at negative voltages. 
For the purposes of this study, it is useful to describe pore 
block in terms of a reduction of average unitary chloride 
current (i) flowing through the open pore (Fig. 1 C; com- 
pare blue and green to red symbols) , as observed in heav- 
ily filtered (10-Hz) records. Plotting fractional unitary 
currents in the presence of NPPB or MOPS - against 
membrane voltage (Fig. 1 D) revealed, despite different 
affinities, similar voltage dependence of pore block: ef- 
fective valences (from Boltzmann fits; Fig. 1 D, solid lines) 
were -0.51 ± 0.03 for NPPB and -0.45 ± 0.03 for MOPS - 
(compare to Ishihara and Welsh, 1997), suggesting that 
both blocking sites sense ^50% of the transmembrane 
electrical field. 

A convenient macroscopic pore block assay is af- 
forded by the nonhydrolytic El 37 IS CFTR mutant, 
which lacks the catalytic glutamate in site 2. These chan- 
nels open rapidly in response to ATP but stay open for 
tens of seconds (Vergani et al., 2003), consistent with 
defective ATP hydrolysis. Because they are active in rest- 
ing cells (compare to Zhou et al., 2010), excision of 
E1371S multichannel patches into an ATP-free solution 
typically uncovers large, slowly decaying currents de- 
void of gating fluctuations (Fig. 1, E, G, I, and K). As 
they flow through virtually locked open channels (P G of 
^1), responses of such currents to brief (2-5-s) applica- 
tions of blockers reflect pure pore block. Indeed, brief 
applications of 210 pM NPPB (Fig. 1 E) or 20 mM 
MOPS - (Fig. 1 G) at various voltages yielded voltage- 
dependent current block (Fig. 1 , F and H, closed sym- 
bols and lines) , which essentially replicated the effects 
on unitary current amplitudes (Fig. 1, F and H, open 
symbols and dotted lines; replotted from Fig. 1 D). 
Using this convenient methodology, we measured pore 
block at varying [NPPB] (Fig. 1 1) and [MOPS ] (Fig. 1 K) 
at a fixed voltage of — 120 mV, constructed dose-response 
curves (Fig. 1,J and L), and obtained apparent Revalues 
of ^20 pM for NPPB (compare to Zhang et al., 2000) 



and ^8 mM for MOPS (compare to Ishihara and 
Welsh, 1997). 

Comparison of effects on macroscopic and unitary 
currents reveals potent gating stimulation by NPPB but 
not MOPS 

We next examined the effects of NPPB and MOPS - on 
macroscopic currents of WT CFTR channels opening 
and closing at steady state in saturating (2 mM) ATP. At 
— 120 mV, the application of increasing [NPPB] (Fig. 2 A) 
or [MOPS ] (Fig. 2 C) inhibited these currents in a 
dose-dependent manner, and fractional inhibition by 
MOPS - (Fig. 2 D, closed symbols; of ^15 mM) was 
reasonably explained by its effect on unitary currents 
(Fig. 2 D, open symbols; replotted from Fig. 1 L). In 
contrast, for NPPB (Fig. 2 B), at any given concentra- 
tion, the fractional reduction of macroscopic WT 
steady-state currents (closed symbols) was milder than 
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Figure 2. Discrepancy between macroscopic and unitary effects 
on WT CFTR quantifies gating stimulation. (A and C) Dose- 
dependent inhibition by cytosolic NPPB (A) and MOPS - (C) 
of steady-state macroscopic WT CFTR currents in 2 mM ATP at 
— 120 mV. (B and D) Dose dependence of fractional currents 
in (B) NPPB and (D) MOPS" at -120 mV (closed symbols), 
and Hill fits (solid lines); open symbols and dotted lines show 
dose dependence of pore block, replotted from Fig. 1 (J and L). 
(E) Gating stimulation by NPPB and MOPS - at -120 mV; 
Po/Po; conuoi was calculated as the ratio (///control) '-{i/ ^control J ■ 
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the fractional reduction of unitary currents (open sym- 
bols; replotted from Fig. 1 J) , yielding an apparent Kj of 
^85 uM (compare to Zhang et al., 2000). 

The macroscopic current (/) is given by /= N ■ i ■ P„, 
where N is the number of channels in the patch, i is 
understood as average unitary current during a burst, 
and P Q is the fraction of time a channel spends in the 
bursting state (bursting probability) . Therefore, a dis- 
crepancy between fractional effects of a drug on i and / 
(Fig. 2 A, red arrow) indicates that the drug must also 
simultaneously affect P„. The fractional effect on P Q 
(Fig. 2 E) is given by the ratio of normalized macro- 
scopic and unitary currents. A marginal increase in P Q 
(^50%) was observed in very high (>20 mM) MOPS" 
(Fig. 2 E, green symbols) . But a much more potent 
gating stimulation was seen with NPPB, at submilli- 
molar concentrations, with P„ increased up to fourfold 
(Fig. 2 E, blue symbols): a close-to-maximal stimulation, 
considering that P„ is M).2 under control conditions 
(see Fig. 5 B). 



NPPB speeds opening and slows closure of WT 
CFTR channels 

To understand the mechanism by which NPPB increases 
P 0 , we systematically studied its effects on the rates of 
individual steps in the gating cycle (Fig. 3 C, cartoon) . 
We first examined effects on the rate of normal, hydro- 
lytic closure of WT CFTR (Fig. 3 C, cartoon, purple 
arrow; r oc ) by comparing macroscopic current relax- 
ation rates (at —120 mV) upon sudden ATP removal in 
the absence or presence of 210 uM NPPB (Fig. 3 A). 
Single-exponential fits (Fig. 3 A, solid lines, time con- 
stants indicated) revealed approximately fourfold slower 
closing rates in the presence of 210 pM NPPB relative to 
control (Fig. 3 C, blue vs. red bar) . In contrast, in simi- 
lar experiments, 20 mM MOPS" (Fig. 3 B) failed to alter 
channel closing rate (Fig. 3 C, green bar). 

To evaluate a potential effect of NPPB on opening 
rate, we exposed WT CFTR channels gating at steady 
state in 2 mM ATP to 210 pM NPPB (Fig. 3 D) . Although 
in similar experiments MOPS - application and removal 
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Figure 3. Effects of NPPB and MOPS" on WT CFTR opening and closing rate. (A and B) Macroscopic WT CFTR currents at - 120 mV 
elicited by brief exposures to 2 mM ATP in the absence or presence of (A) 210 uM NPPB or (B) 20 mM MOPS - ; colored lines, single- 
exponential fits (t, time constants). (C) Macroscopic closing rates (bars, r 0 c = 1/t) in the absence (red) and presence of NPPB (blue) or 
MOPS - (green) . (D and E) Macroscopic WT CFTR currents in 2 mM ATP at - 120 mV, and brief exposure to (D) 210 uM NPPB or (E) 
20 mM MOPS - . Colored lines in D, single-exponential fits. (F) Sums of opening and closing rates in the absence and presence of NPPB 
(1/t from D); opening rates (r C o, bottom vertical arrows) are estimated by subtraction of closing rates (r 0 c, top vertical arrows; see C) 
from the sums. Cartoons in C and F (also in Figs. 4, B and E, and 5, B and E) depict simplified cyclic CFTR gating model: cyan, TMDs; 
green, NBD1; blue, NBD2; yellow, ATP; red, AD P. Purple arrows highlight the pathway under study. 
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caused simple monophasic current relaxations (Fig. 3 E) , 
both the addition and removal of NPPB elicited tripha- 
sic responses (Fig. 3 D) , attesting to the dual effects of 
this compound. Upon NPPB application, instantaneous 
(t of ^20 ms, reflecting solution exchange time) pore 
block was followed by a partial current recovery, as a 
larger pool of channels was drawn into the open burst 
state. Upon NPPB withdrawal, instantaneous relief from 
pore block revealed this larger pool of open channels in 
the form of a large current overshoot, which then re- 
laxed back to the pre-application steady-state current 
level. The rate of macroscopic current relaxation after a 
sudden change in gating parameters reflects the sum of 
the average microscopic opening and closing rates in 
the new conditions. Interestingly, single-exponential fits 
to the current time courses after block (Fig. 3 D, blue fit 
line) and unblock (red fit line) yielded similar time 
constants, suggesting that the sum of opening and clos- 
ing rates in the presence (Fig. 3 F, blue bar) and ab- 
sence (Fig. 3 F, red bar) of 210 pM NPPB is not very 
different. However, in the absence of NPPB, this sum is 



largely dominated by the closing rate (compare Fig. 3 F, 
red bar, to Fig. 3 C, red bar). In contrast, in the pres- 
ence of NPPB, closing rate is slowed to ^1 s _1 (Fig. 3 C, 
blue bar) ; thus, the sum (Fig. 3 F, blue bar) must be dom- 
inated by an increased opening rate (see also Fig. 6 D) . 
These data suggest that NPPB increases channel open- 
ing rate, which, under our conditions of saturating ATP, 
reflects the rate of step Ci— >Oj (Fig. 3 F, cartoon, purple 
arrow; r C o) • 

In a nonhydrolytic CFTR mutant, NPPB speeds gating 
but does not affect P c 

For most ion channels, potentiators stabilize open states. 
If NPPB acted by stabilizing state Oi relative to Ci, 
then, in addition to speeding opening (see Fig. 3 F), it 
might be expected to slow the reverse of the opening 
step, i.e., nonhydrolytic closure (Fig. 4 B, cartoon, pur- 
ple arrow; k\) . To test this, we studied the closing rate of 
K1250A CFTR channels, in which mutation of the 
NBD2 Walker A lysine abrogates ATP hydrolysis at site 2 
(Ramjeesingh et al., 1999) and reduces gating to reversible 




Figure 4. Effects of NPPB and MOPS - on gating rates under nonhydrolytic conditions. (A) Macroscopic K1250A CFTR current at 
— 120 mV elicited by exposures to 10 mM ATP in the absence or presence of blockers. Colored lines, single-exponential fits (t, time 
constants). (B) Macroscopic closing rates (bars, 1/t) in the absence (red) and presence of NPPB (blue) or MOPS - (green) quantify 
effects on rate k\ (cartoon, purple arrow). The K1250A mutation (B and E, cartoons, red stars) disrupts ATP hydrolysis in site 2 (red 
cross). (C) Macroscopic K1250A CFTR current elicited by 10 mM ATP at -40 mV, prolonged exposure to 210 uM NPPB of channels 
gating at steady state, and brief exposure to NPPB of surviving locked-open channels after ATP removal ( 10-s yellow box, expanded in 
inset). Bracketing brief applications of NPPB to small residual CFTR currents were used to estimate zero-current level (dotted line). 
(D) Fractional K1250A CFTR currents at —40 mV in 210 pM NPPB applied during steady-state gating (red bar) or in the locked-open 
state (yellow bar). (E) Effect of NPPB on the closed-open equilibrium (cartoon, purple double arrow). Fractional effect on P 0 for 
K1250A CFTR (blue bar) was calculated as in Fig. 2 E. 
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Ci<-»Oi transitions (Fig. 4 B, cartoon). Consistent with 
previous reports (Gunderson and Kopito, 1995; Vergani 
et al., 2003; Csanady et al., 2010), macroscopic closing 
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Figure 5. Effects of NPPB and MOPS" on WT CFTR micro- 
scopic steady-state gating parameters. (A) Currents from single 
WT CFTR channels at - 120 mV in 2 mM ATP ± blockers; band- 
width, 100 Hz. (B) Mean burst (T b ) and interburst (T ib ) durations 
(bars) in 2 mM ATP (red), 2 mM ATP plus 20 uM NPPB (blue), 
and 2 mM ATP plus 20 mM MOPS" (green); closing rate (car- 
toon, purple arrow) is 1 / T b . Bars were generated from data origi- 
nating from 15, 10, and 9 patches, respectively, for the control, 
NPPB, and MOPS conditions; error bars represent SEM. (C and D) 
Dwell-time histograms of burst durations of WT CFTR chan- 
nels at -120 mV in 2 mM ATP (C; 2,999 bursts, pooled from 15 
patches) or 2 mM ATP plus 20 uM NPPB (D; 2,069 bursts, pooled 
from 10 patches) and maximum likelihood fits (solid lines) to 
the scheme cartooned in E (with ki fixed to zero). (E) Estimates 
of rates k x (step Oi->0 2 ) and fe (step 0»->C 2 ) in ATP and ATP 
plus 20 uM NPPB. Asymmetric error bars represent 0.5-unit log- 
likelihood intervals. 



rate of K1250A upon removal of bath ATP (Fig. 4 A) was 
^100 times slower than for WT channels (Fig. 4, A, red 
fit line, and B, red bar; compare to Fig. 3, A and C) . 
Whereas 20 mM MOPS" did not affect nonhydrolytic 
closure (Fig. 4, A, green fit line, and B, green bar) , 210 uM 
NPPB unexpectedly accelerated it by two- to three- 
fold (Fig. 4, A, blue fit line, and B, blue bar) . As an alter- 
native approach for studying effects on nonhydrolytic 
closure, we also compared unlocking rates in the ab- 
sence and presence of NPPB for WT CFTR channels 
locked in the open state by exposure to a mixture of 
ATP plus pyrophosphate (Gunderson and Kopito, 1994; 
Tsai et al., 2009). Consistent with its effect on K1250A 
closing rate (Fig. 4, A and B), 210 uM NPPB also accel- 
erated the unlocking rate of pyrophosphate-locked WT 
CFTR channels by two- to threefold (Fig. S2). 

Because NPPB accelerated both forward (Fig. 3 F) 
and backward (Fig. 4 B) transitions of the Ci<->Oi step, 
we examined whether it affects the equilibrium con- 
stant between those two states (Fig. 4 E, cartoon, pur- 
ple double arrow; K<. q ), i.e., P„ for K1250A channels. 
Even prolonged exposure to 210 uM NPPB of K1250A 
channels gating at steady state in 10 mM ATP (Fig. 4 C; 
note V m = — 40 mV) failed to elicit biphasic responses 
like those seen for WT channels (see Fig. 3D). More- 
over, fractional current inhibition at steady state was 
identical to that instantaneously observed upon brief 
application of NPPB long after ATP removal (Fig. 4 C, 
yellow box, expanded in inset) . The latter maneuver 
measures pure pore block of surviving open channels, 
i.e., fractional reduction of i (see Fig. 1, E, G, I, and K). 
The very similar fractional NPPB effects under those 
two conditions (Fig. 4 D, red vs. yellow bar; compare to 
Fig. 2 B) revealed no fractional change in P„ (Fig. 4 E, 
blue bar) . Thus, in contrast to its effect on hydrolytic 
gating (Fig. 2 E), NPPB does not affect P Q for steady- 
state nonhydrolytic gating. Given the observed in- 
crease in nonhydrolytic closing rate (Fig. 4, A and B), 
this result indicates that NPPB must also speed the 
opening rate of K1250A CFTR, just as it does for 
WT (Fig. 3 F) . The magnitude of the NPPB effects 
on opening and closing rates must be similar, such 
that the Ci<->Oi equilibrium constant is not affected; 
i.e., by reducing the height of the energetic barrier 
separating the two states, NPPB acts as a catalyst for 
this step. 

NPPB slows WT CFTR closure by slowing the ATP- 
hydrolysis step 

NPPB slowed the hydrolytic closing rate of WT CFTR 
(Fig. 3, A and C) , which reflects sequential transition 
through the prehydrolytic Oi and post-hydrolytic O2 
states (Fig. 3 C, cartoon, purple arrow; r 0 c) • Dissection 
of the rates of the two sequential steps requires maxi- 
mum likelihood fitting of the distributions of open 
burst durations (Csanady et al., 2010). To achieve this, 
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we analyzed steady-state gating of single WT CFTR chan- 
nels under control conditions, in 20 uM NPPB (unitary 
currents in 210 uM NPPB could not be resolved at 
-120 mV [Fig. 1 A], whereas 20 uM NPPB did not com- 
promise reliable event detection) , or in 20 mM MOPS - 
(Fig. 5 A). 20 mM MOPS - did not significandy affect 
either mean burst (308 ± 60 ms; n = 9) or mean inter- 
burst durations (Fig. 5 B, green bars) . However NPPB, 
even at this low concentration, significantly prolonged 
mean burst duration (from 277 ± 24 ms [n = 37] to 412 ± 
56 ms [n = 15]; P = 0.013), i.e., slowed closure; reduction 
of mean interburst duration was not significant (Fig. 5 B, 
blue vs. red bars) . 

To determine which of the two sequential closing 
steps is slowed by NPPB, we reconstructed the distribu- 
tions of burst durations for single WT CFTR channels 
gating at — 120 mV in the absence (Fig. 5 C) or presence 
of 20 uM NPPB (Fig. 5 D). Both distributions were 
peaked and fit (Fig. 5, C and D, solid lines) significandy 
better (P = 4 x 10~ 12 and 9 x 10~ 7 , respectively, by the 
log-likelihood ratio test) by the two-parameter (fa, hi) 
irreversible sequential closing mechanism (Fig. 5 E, car- 
toon, purple arrows) than by a single exponential. The 
fit parameters (Fig. 5 E, bars) attested to a significant 
(P < 10~ 4 ) reduction by NPPB of the slow rate fa, which 
represents the rate of step Oi— >02 associated with ATP 
hydrolysis in site 2 (Csanady et al., 2010) . 

Voltage-independent gating stimulation by NPPB suggests 
distinct binding sites for pore block and gating effects 
NPPB was first identified as a potentiator at positive volt- 
ages (Wang et al., 2005) because at negative voltages 
gating stimulation is masked by the pore block, which is 
voltage dependent. To assess whether any of its effects 
on gating are also voltage dependent, we systematically 
studied gating in NPPB at +60 mV and compared NPPB 
effects at +60 mV with those at -120 mV (Figs. 1-4). Be- 
cause at positive voltages unitary current in 210 uM 
NPPB remains ^67% of control (as opposed to M)% at 
— 120 mV; Fig. 1 D), gating stimulation prevails over 
pore block, causing dose-dependent overall stimulation 
of macroscopic WT CFTR current (Fig. 6 A) . Indeed, 
dividing fractional macroscopic currents (Fig. 6 B, open 
squares) by fractional unitary currents (Fig. 6 B, open 
diamonds) revealed a powerful increase in P„ (Fig. 6 B, 
closed circles); the approximately threefold stimulation 
by 210 uM NPPB approached that observed at - 120 mV 
(Fig. 2 E, blue circles). Furthermore, 210 pM NPPB 
slowed closing rate by approximately threefold (Fig. 6 C; 
compare to Fig. 3, A and C) , and similar macroscopic 
relaxation rates upon the addition and removal of NPPB 
attested to a simultaneous increase in opening rate 
in the presence of NPPB (Fig. 6 D; compare to Fig. 3, 
D and F). Finally, at +60 mV, 210 pM NPPB increased 
the nonhydrolytic closing rate of K1250A CFTR by ap- 
proximately threefold (Fig. 6 E) Just as it did at — 120 mV 




Figure 6. Gating stimulation by NPPB is largely voltage inde- 
pendent. (A) Dose-dependent stimulation by NPPB of steady- 
state macroscopic WT CFTR currents in 2 mM ATP at +60 mV. 
(B) Dose-response curves for fractional effects of NPPB on steady- 
state macroscopic current I (open squares; obtained from A), on 
unitary current i (open diamonds; obtained as in Fig. II), and on 
P 0 (closed circles; obtained as in Fig. 2 E). (C; left) Macroscopic 
WT CFTR current at +60 mV elicited by brief exposures to 2 mM 
ATP in the absence or presence of 210 uM NPPB; colored lines, 
single-exponential fits (t, time constants). (Right) Macroscopic 
WT closing rates (bars, 1/t) in the absence (red) and presence 
of NPPB (blue). (D; left) Macroscopic WT CFTR current in 
2 mM ATP at +60 mV and brief exposure to 210 uM NPPB; colored 
lines, single-exponential fits (t, time constants). (Right) Sums of 
opening and closing rates (bars, 1/t) in the absence and pres- 
ence of NPPB; opening rates (r C o, bottom vertical arrows) are es- 
timated by subtraction of closing rates (r 0 c, top vertical arrows) 
from the sums. (E; left) Macroscopic K1250A CFTR current at 
+60 mV elicited by exposures to 10 mM ATP in the absence or 
presence of 210 pM NPPB; colored lines, single-exponential fits 
(t, time constants). (Right) Macroscopic K1250A closing rates 
(bars, 1/t) in the absence (red) and presence (blue) of NPPB. 
(F; top) Currents from a single WT CFTR channel gating at 
+60 mV in 2 mM ATP before (top trace) and during (bottom trace) 
exposure to 210 pM NPPB; bandwidth, 100 Hz. (Bottom) Mean 
burst (left) and interburst (right) durations (bars) of WT CFTR 
in 2 mM ATP (red) and 2 mM ATP plus 210 pM NPPB (blue) 
at +60 mV. Bars were generated from 1 6 segments of record in 
ATP and 15 segments in ATP plus NPPB originating from nine 
patches; error bars represent SEM. 
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(Fig. 4, A and B) . Because at +60 mV unitary currents 
remained resolvable even in the presence of 210 pM NPPB 
(Fig. 6 F) , the conclusions drawn from nonstationary 
macroscopic analysis regarding effects on WT CFTR gat- 
ing (Fig. 6, C and D) could be corroborated in steady- 
state single-channel recordings: as expected, 210 uM 
NPPB prolonged mean burst duration but shortened 
mean interburst duration, both by approximately three- 
fold (Fig. 6 F, blue bars) . 

Stimulation of P Q by NPPB is enhanced for poorly active 
AF508 or unphosphorylated WT CFTR channels 
In contrast to the three- to fourfold stimulation by 
NPPB of P 0 for WT CFTR observed here, Wang et al. 
(2005) reported 10-15-fold stimulation of AF508 CFTR 
by NPPB. Indeed, we found that gating stimulation of 
low temperature-rescued AF508 CFTR by 210 pM 
NPPB is so robust that it overrides the pore block effect 
at all voltages, yielding overall current stimulation even 
at —120 mV (Fig. 7, A and C, yellow bar), and greater 
than fivefold stimulation at +60 mV (Fig. 7, B and C, 



Figure 7. Stimulation of AF508 CFTR gating by NPPB. 
(A and B) Quasi-macroscopic currents of low temperature- 
rescued, prephosphorylated AF508 CFTR channels in 
2 mM ATP at - 1 20 mV (A) and +60 mV (B) , and brief expo- 
sures to 210 uM NPPB. Note overall current stimulation by 
NPPB, and extreme current overshoot upon its removal, at 
-120 mV. (C and D) Fractional effects of 210 uM NPPB on 
(C) steady-state macroscopic currents and (D) P„ of AF508 
CFTR at — 120 mV (yellow bars) and +60 mV (cyan bars) . 

cyan bar) . Accounting for the pore block effect, these 
results suggest ^20- and ^8-fold stimulation of P Q at 
— 120 and +60 mV, respectively (Fig. 7 D). 

Interestingly, Wang et al. (2005) also reported a de- 
pendence of NPPB stimulation on the phosphorylation 
state of WT CFTR. Whereas all the data in Figs. 2-6 
were obtained in the relatively stable, partially de- 
phosphorylated state of CFTR, after removal of PKA 
(Csanady et al., 2010), we indeed observed an ^10-fold 
stimulation of the small current, carried by low P., activ- 
ity of unphosphorylated WT CFTR channels, when 210 
pM NPPB was applied at +60 mV before exposure of the 
patch to PKA (Fig. 8, A, first application of NPPB [ex- 
panded in inset], and B, left). In contrast, in the same 
patch, 210 pM NPPB had little effect on overall current 
at +60 mV when applied in the presence of PKA, i.e., to 
fully phosphorylated channels that are gating at high P Q 
(Fig. 8, A, second application of NPPB, and B, right) , 
whereas soon after PKA removal, it increased by approxi- 
mately twofold the current carried by partially dephos- 
phorylated channels, which are gating at intermediate 
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Figure 8. Phosphorylation dependence of NPPB stimu- 
lation of WT CFTR reflects nonlinear dependence of P 0 
on opening and closing rates. (A) Macroscopic WT CFTR 
currents in 2 mM ATP at +60 mV, and applications of 
210 pM NPPB before exposure to (current scale magnified in 
inset) , in the presence of, and after removal of 300 nM PKA 
catalytic subunit. (B and C) Fractional effects of 210 pM 
NPPB on (B) steady-state macroscopic currents and (C) 
P 0 of unphosphorylated (pre-PKA), fully phosphorylated 
(in PKA), and partially dephosphorylated (post-PKA) 
WT CFTR channels at +60 mV. (D) Table summarizing ex- 
pected changes in gating parameters assuming a constant 
(phosphorylation-independent) fourfold increase in T b and 
fourfold decrease in T ib in the presence of 210 pM NPPB, 
as measured in this study for the post-PKA condition at 
— 120 mV (Fig. 3; slightly smaller, approximately threefold, 
effects were measured at +60 mV; Fig. 6) . Control parameters 
(left column) reflect typical values measured for WT CFTR in 
this and previous studies (e.g., Csanady et al., 2005; Szollosi 
et al., 2010); the post-PKA values in NPPB (right column, 
bottom row) reflect the values measured in Fig. 3. 
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Figure 9. Anionic form of NPPB is responsible for the observed 
gating effects. (A) Titration curve of NPPB; 10 ml 100 uM NPPB 
(in H 2 0) was titrated with 10-pl aliquots of either 10 mM NaOH 
(closed symbols; n = 3) or 10 mM HQ (open symbols; n = 2, plot- 
ted on negative side of abscissa) . Solid line is a fit to the titration 
curve of a monoprotic weak acid, corrected for the presence of 
10 mM C0 2 (see Materials and methods); fitted pK,, is plotted. 
(B) Structures of protonated (zwitterionic) and deprotonated 
(anionic) forms of NPPB; the proton released with a pK„ of ~5.6 
is highlighted in red; dotted line depicts suggested hydrogen 
bond. (C) Segments of unitary current at +60 mV from a single 
locked-open E1371S CFTR channel in symmetrical ^140-mM 
Cl~ at three different cytosolic pH values, with and without cyto- 
solically applied 210 pM NPPB, displayed at two bandwidths 
(f c , filter corner frequency) . Horizontal lines depict closed current 
level and average current through a bursting channel. (D) Ap- 
parent unitary amplitudes (of heavily filtered currents) , normal- 
ized to those obtained in the absence of NPPB at the respective 
pH (bars) . Table lists concentrations of anionic and uncharged 
NPPB for the tested conditions, calculated using a pK,, of 5.6. 
(E and F) Stimulation at +60 mV of steady-state macroscopic WT 
CFTR currents in 2 mM ATP (red bars) by the application of 



P„ (Fig. 8, A, last NPPB application, and B, middle) . Ac- 
counting for the pore block, these effects of NPPB on 
overall WT CFTR current at +60 mV correspond to 
~15-, ^3-, and <2-fold stimulation of P Q of poorly, par- 
tially, and fully phosphorylated CFTR (Fig. 8 C) . These 
data confirm the findings of Wang etal. (2005) and sug- 
gest that NPPB stimulates P (> much more efficiently 
when the starting P Q value is low. 

Anionic form of NPPB is responsible for the observed 
gating stimulation 

At our bath pH of 7.1, the anionic, deprotonated form of 
NPPB dominates, and only a small fraction of the com- 
pound is protonated, uncharged. Although the pore block 
is no doubt attributable to the anionic form of NPPB 
(NPPB ), we wanted to identify whether it is this anionic 
form or perhaps the less prevalent protonated form 
(NPPB-H) that is responsible for gating stimulation. 

Because we found no reliable data in the literature for 
the pK„ value of NPPB (most studies assumed a pK a of 
^4.5 by referring either to the study of Wangemann et 
al., 1986, in which the pK,, was not measured, or to that of 
Walsh and Wang, 1998, in which no data were shown), 
we determined the latter by electrometric titration. 
Because of the limited solubility of NPPB in water (espe- 
cially at low pH), titration by NaOH (Fig. 9 A, closed 
circles) and HC1 (Fig. 9 A, open circles) could be per- 
formed only for a dilute (100-uM) aqueous solution of 
the compound, restricting reliable determination of po- 
tential pK,, values to the range of 4-10. (At this low 
concentration, a weak acid with a pK„ of <4 is largely 
dissociated, whereas a site with a pK„ of >10 remains 
largely protonated even after the addition of 1 mol/ mol 
NaOH.) The data were well fit (Fig. 1 A, solid line) by a 
titration curve expected for a monoprotic weak acid and 
revealed a single pK,, value of ^5.6 in the range of 4—10. 

To deconvolve gating and permeation effects of NPPB 
at various pH, we first characterized pore block at 
+60 mV by 210 pM of total [NPPB] applied at bath pH 
values of 8.1, 7.1, and 6.1 (Fig. 9 C). Whereas average 
unitary currents, determined from heavily filtered cur- 
rent traces, decreased to ^67% of control when 210 pM 
NPPB was applied at pH 7.1 or 8.1 (Fig. 9 D, left; com- 
pare to Fig. 1 D), at pH 6.1 the pore block was signif- 
icantly (P = 8 x 10~ 4 ) milder, yielding a fractional 
decrease in unitary current of only ~15% (Fig. 9 D, 
third bar) . This is consistent with a significantly de- 
creased concentration of NPPB - at this low pH, as 



(E and F) 210 uM NPPB (blue bars) at bath pH values of 6.1, 7.1, 
and 8.1 (gray bars), or of (F) 21 pM NPPB at a bath pH value of 
6.1. (G) pH dependence of fractional effects of NPPB on steady- 
state macroscopic current / (blue bars; obtained from E and F), 
and on P„ (dark blue bars; P 0 /P o; control was calculated as the ratio 
{I/Icamo\)-{i/ Control))- (Table) Concentrations of anionic and un- 
charged NPPB for the tested conditions, assuming pK„ = 5.6. 
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predicted by the measured pK a value of 5.6 (see table 
in Fig. 9D). 

We next assayed macroscopic current stimulation at 
+60 mV by 210 pM of total NPPB at various pH (Fig. 9 E) . 
Fractional current stimulation at pH 6.1 was not en- 
hanced relative to that measured at our control pH of 
7.1 (Fig. 9 E and middle blue bars in G), despite an 
^10-fold higher concentration of uncharged NPPB-H 
at the lower pH (Fig. 9 G, table); in fact, the calculated 
effect on P c was significantly blunted at pH 6.1 (Fig. 9 G, 
middle dark blue bars). Furthermore, no significant 
stimulation of P„ was observed when applying 21 pM 
of total NPPB at pH 6.1 (Fig. 9 F and right side of G), 
a condition in which [NPPB ] is selectively lowered, 
whereas [NPPB-H] remains comparable to that calcu- 
lated for 210 pM of total NPPB at pH 7.1 (Fig. 9 G, 
table) . These results rule out uncharged NPPB-H as the 
active form responsible for gating stimulation and sug- 
gest that the latter is largely caused by anionic NPPB . 

In principle, the reduced fractional stimulation of P Q 
by 210 pM of total NPPB at pH 6.1 could be taken as an 
indication that uncharged NPPB-H is entirely inactive; 
however, we refrain ourselves from such a conclusion, 
as fractional current stimulation was also blunted at pH 
8.1 (Fig. 9 E and left blue bar in G) relative to pH 7.1, 
despite slightly higher NPPB - at the higher pH (Fig. 9 G, 
table). The similarly reduced fractional effect on P Q at 
both pH 6.1 and 8.1 (Fig. 9 G, compare first and third 
dark blue bars) suggests that changes in surface proper- 
ties and/ or gating kinetics of the CFTR protein at both 
acidic and basic pH should also be considered (compare 
to Chen etal., 2009). 

DISCUSSION 

Until very recently, treatment of CF has been exclusively 
symptomatic. This has changed with the approval of VX- 
770 (ivacaftor; Vertex Pharmaceuticals; Van Goor et al., 
2009; Ramsey et al., 201 1 ) , the first drug shown to act by 
binding specifically to the CFTR protein. However, to 
date it has been shown to benefit only patients carrying 
at least one G551D allele, which constitute <5% of the 
CF population. Effective drug therapy to treat patients 
carrying AF508 alleles is still lacking. At present, there is 
a wide gap between industrial efforts to obtain drugs tar- 
geting CFTR and academic research aimed at under- 
standing CFTR structure and mechanism at a molecular 
level. Although industrial high-throughput screening 
has clearly met some success (Van Goor et al., 2009; 
Ramsey et al., 2011), the wealth of information that has 
emerged from basic research has yet to impact drug dis- 
covery. In this study, we have identified two strategic 
points in CFTR's unique functional cycle that are emi- 
nently suited as intervention points for altering CFTR 
activity: the channel opening step, and the step that rate 
limits channel closure (hydrolytic Oi— >0% step). Our 



dissection of the mechanism of action of NPPB, pre- 
sented here, shows how affecting these transitions can 
powerfully influence P G of WT and mutant (including 
AF508) CFTR channels, and offers a route to rational 
design of drugs targeting CFTR. 

Our interpretation builds on the nonequilibrium 
nature of the CFTR gating cycle (see the simplified cy- 
clic gating model depicted throughout the insets of 
Figs. 3-5) , which is rooted in its evolutionary descent 
from an ancestral ABC transporter. First proposed by 
Vergani et al. (2003), this model was later tested and 
refined (Vergani et al., 2005; Csanady et al., 2006, 2010; 
Bompadre et al., 2007; Tsai et al., 2010) before becom- 
ing generally accepted (Gadsby et al., 2006; Chen and 
Hwang, 2008; Hwang and Sheppard, 2009; but compare 
to Aleksandrov et al., 2009). Recently, various exten- 
sions of this scheme have been suggested to explain the 
rare openings (P Q of ^0.004) , in the absence of ATP, of 
unliganded WT (Szollosi et al., 2010) or G551D CFTR 
(Bompadre et al., 2007), enhancement of such activity 
by cytosolic loop mutations (Wang et al., 2010), and 
[ATP] -dependent prolongation of open burst dura- 
tions in mutants (Jih et al., 2012b) or by drugs (Jih and 
Hwang, 2013). Although these extended models differ 
in detail (Szollosi et al., 2010; Kirk and Wang, 2011; Jih 
et al., 2012a), they all contain, as a common core, the 
basic cycle depicted in our figure insets, which describes 
the majority of gating events for both WT CFTR and 
many CFTR mutants, including AF508 (Jih et al., 2011; 
but compare to Wang et al., 2005, for ATP-independent 
G551D). We will therefore use the core cyclic model as 
a framework for discussing our results. 

Voltage-dependent pore block and gating stimulation 
are distinct NPPB effects 

Our study clearly distinguishes the effects of NPPB on 
gating from its voltage-dependent pore-blocker effects. 
Although stimulation by NPPB was originally only noted 
at positive potentials (Wang et al., 2005), careful com- 
parison of effects on unitary versus macroscopic cur- 
rents demonstrates robust stimulation of P Q even at 
— 120 mV (Fig. 2) . In fact, because of its dual action on 
both i and P„, fractional reduction of macroscopic cur- 
rent by NPPB is not a good measure of pore block, un- 
less observed on "locked-open" channels, which are not 
gating (P„ of ^T), such as surviving E1371S (Fig. 1, 
E and I) or K1250A (Fig. 4 D, inset) channels after the 
removal of ATP. This might explain previously reported 
lower values (-0.35, Zhang et al., 2000; -0.2, Zhou 
et al., 2010) for the effective valence of the NPPB block 
than found here (approximately —0.5; Fig. 1, D-F). Of 
note, Ai et al. (2004) found an effective valence of ap- 
proximately — 0.5 for CFTR pore block by anthracene-9- 
carboxylate (9-AC), a compound that also stimulates P Q 
of CFTR, suggesting that NPPB, MOPS, and 9-AC likely 
share a common binding site in the pore; possibly, gating 
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stimulation by NPPB and 9-AC might also involve a 
shared binding site. 

We show here that all gating effects of NPPB are in- 
deed largely voltage independent at a microscopic level 
(Fig. 6, A-F; except for slighdy smaller effects at +60 mV) , 
as expected if gating and pore block effects involve 
distinct binding sites. Accordingly, although the appar- 
ent of NPPB for pore block was ^20 uM at - 1 20 mV 
(Fig. 1 J; replotted in Fig. 10, open diamonds), the ap- 
parent K1/2 for slowing hydrolytic closure of WT CFTR 
(Fig. 10, closed circles) was at least four times higher 
(M)0 uM) , and a tentative fit to the dose-response curve 
for acceleration of nonhydrolytic closure (measured for 
the K1250A mutant; Fig. 10, closed squares) suggested a 
similar K 1/2 of ^100 uM (although reliability of the lat- 
ter fit is limited by the uncertainty of its asymptotic 
value). Consistent with our interpretation of distinct 
binding sites for pore block and gating effects, replace- 
ment of the negatively charged carboxyl moiety of 
NPPB by an uncharged amide group eliminated pore 
block and yielded a pure potentiator compound, which 
stimulated gating of both WT and AF508 CFTR at all 
voltages (Wang et al., 2005). 

Anionic form of NPPB stimulates gating 
Although the uncharged amide analogue of NPPB stim- 
ulates CFTR (Wang et al., 2005), the gating effect of 
NPPB was not enhanced at acidic pH, ruling out un- 
charged, protonated NPPB-H, which is present at low 
micromolar concentrations at our standard bath pH of 
7.1, as the cause of the observed gating stimulation. We 
therefore conclude that the stimulatory effect observed 
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Figure 10. Dose dependence of NPPB gating effects. Macro- 
scopic closing rates of WT (closed circles) and K1250A (closed 
squares) CFTR in the presence of various cytosolic [NPPB] , mea- 
sured at — 120 mV using the protocols shown in Figs. 3 A and 4 A, 
respectively; leftmost data points represent control values in the 
absence of NPPB. The 50- and 100-uM data points for K1250A 
were measured at —40 mV. Solid lines are fits to the equation 
r oc ([NPPB]) = r 0 + (r«, - r„) ([NPPB]/ ( [NPPB] + K 1/2 )), with r„ 
fixed to zero for WT but left free for K1250A; K l/2 values are plot- 
ted. Open diamonds and dotted line depict fractional pore block 
at —120 mV, replotted from Fig. 1 J. 



in this study is largely caused by anionic NPPB . Of 
note, NPPB-H is not expected to be a good mimic of the 
amide analogue. This is because the single pK„ value of 
^5.6, which we could detect in the range of pH 4—10 
(Fig. 9 A), is unlikely to reflect titration of the carboxyl- 
ate group, but rather that of the secondary amino group 
of NPPB (Fig. 9 B). Indeed, in 3-nitrobenzoic acid, the 
pK a of the carboxylate is ^3.5, and for 2-aminobenzoic 
acid, the two pK„ values are ^2.1 (for the carboxylate) 
and ^5.0 (for the amino group). (Referenced pK^, val- 
ues were taken from a table available at http://www. 
zirchrom.com/organic.htm.) The unusually low pK„ 
for the amino group in the latter compound is caused 
by the electron-withdrawing effect of the aromatic ring 
(compare, pK„ is ^4.7 for aniline and ^5.1 for the sec- 
ondary amino group of JV-ethyl aniline), whereas the 
unusually low pK,, of its carboxylate group is partly 
caused by a hydrogen bond between the adjacent amino 
group (donor) and a carboxylate oxygen (acceptor); 
such a hydrogen bond is likely to form also in NPPB 
(Fig. 9 B, dotted line). Thus, we interpret the measured 
pK„ value of ^5.6 for NPPB to reflect that of the amino 
group, whereas that of the carboxylate must be too low 
(<3) to be detected in a dilute (100-uM) solution. This 
assignment predicts that NPPB-H, the uncharged form 
of NPPB that accumulates at pH 6.1 (Fig. 9 G, table) , is 
in fact a zwitterion, unlike its amide analogue. 

One explanation for the smaller stimulation of P D by 
NPPB at low pH would be if zwitterionic NPPB-H did 
not act as a stimulator at all. However, because the ef- 
ficiency of gating stimulation declined both at acidic 
and basic pH (Fig. 9 G, dark blue bars) , it is likely that 
changes in the surface charge distribution of the CFTR 
protein, as well as its gating kinetics, are at least partly 
responsible for the decreased responses at extreme pH. 
Indeed, in experiments performed at 37°C, Chen et al. 
(2009) found that P 0 of WT CFTR decreases monotoni- 
cally with increasing pH, within the range of 6.3-8.3. 
The approximately twofold higher P 0 at pH 6.3 com- 
pared with 8.3 was largely caused by approximately 
three- to fourfold longer mean open burst durations at 
the lower pH. Although at 25°C steady-state values 
of macroscopic currents appeared less sensitive to pH 
(Fig. 9 E) , current relaxation rates upon the addition of 
ATP, which reflect the sum of the opening and closing 
rate, were indeed visibly slower at pH 6.1, but faster at 
pH 8.1, when compared with those observed at pH 7.1 
(compare to Fig. 9 E) , reporting profoundly altered gat- 
ing kinetics at both extreme values of pH. In addition to 
the above instantaneous effects of pH on gating kinet- 
ics, we also typically observed a progressive decline in 
macroscopic CFTR currents upon prolonged exposure 
to pH 8.1 (see Fig. 9 E), suggesting that the CFTR pro- 
tein is either dephosphorylated more rapidly, or other- 
wise progressively inactivates at very high pH. These 
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limitations impede us from concluding whether zwitter- 
ionic NPPB-H can act as a gating stimulator or not. 

NPPB catalyzes channel opening 

NPPB stimulates P u by two distinct mechanisms: one 
increasing opening (and nonhydrolytic closing) rate, the 
other slowing the microscopic rate of hydrolysis (Oi— >Oo 
transition). The first action is a stabilization by ^1 AT 
of the Ci<->-Oi transition state, a bona fide catalyst effect, 
which results in comparable (approximately threefold) 
acceleration of "forward" channel opening rate (Figs. 3 F 
and 6, D and F) and "backward" nonhydrolytic closing rate 
(Figs. 4 B and 6 E) . At first glance, it might seem surpris- 
ing that a catalyst should enhance the P„ of an ion chan- 
nel. Indeed, for an equilibrium mechanism (Fig. 11 A, 
left, red rates) , a catalyst, by lowering the free-energy bar- 
rier AG* (Fig. 1 1 B, left, blue vs. red AG profile) , simulta- 
neously speeds both the opening and closing rate by 
the same factor (Fig. 11 A, left, blue rates), causing no 
change in P Q (Fig. 11 C, left). Accordingly, for the non- 
hydrolytic K1250A mutant, which gates at equilibrium, 
NPPB speeds gating transitions but does not alter P Q 
(Fig. 4 E) . In contrast, WT CFTR, because of its unique 
nonequilibrium cycle, enters and exits the open state 
through different pathways, rate-limited by the Q— >Oi 
and the Oi— >Oa transitions, respectively; the rates of these 
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Figure 1 1 . Catalyst effect as a unique strategy to enhance P„ of 
a nonequilibrium ion channel. (A) Simple equilibrium scheme 
(left) and cyclic nonequilibrium scheme (right; compare to car- 
toons in Figs. 3-5), with transition rates in the absence (red) 
and presence (blue) of a catalyst for step G<->0 (left) or Ci<->Oi 
(right). (B) Free energy profiles (not drawn to scale) for the 
schemes in A, in the absence (red) and presence (blue) of the 
same catalyst. (C) Illustration of channel gating patterns, and cal- 
culated P 0 , for the schemes in A, in the absence (red) and pres- 
ence (blue) of the catalyst. 



transitions are determined by the height of the corre- 
sponding energetic barriers (AG*ci _ >oi and AG + oi— >o<i, 
Fig. 11 B, right, red AG profile). Thus, for WT CFTR, 
lowering the Ci<->Oi barrier (Fig. 11 B, right, blue AG 
profile) selectively speeds the opening rate, without 
speeding the closing rate (as long as AG + oi— >a remains 
substantially larger than AG*oi— >oa) • For transition rates 
typical to WT CFTR under our recording conditions 
(Fig. 11 A, right, red rates), this catalyst effect alone of 
NPPB (Fig. 1 1 A, right, blue rates) provides an approxi- 
mate two- to threefold stimulation of P„ (Fig. 1 1 C, right) . 

Thermodynamic studies have provided some insight 
into the nature of the opening transition state, outlin- 
ing it as a high enthalpy structure in which the NBD 
dimer has already formed, but the pore is still closed, 
causing molecular strain at the TMD-NBD interface 
(Csanady et al., 2006; Wang et al., 2010). Recent ABC 
transporter crystals revealed the 3-D structure of this 
interface and how it changes in different functional 
states (Dawson and Locher, 2006; Hohl et al, 2012). For 
CFTR, a homology model, thought to represent a chan- 
nel open state, predicted physical interactions between 
phenylalanine 508 on the NBD1 surface and a cluster of 
aromatic residues in the fourth intracellular TMD loop 
(Serohijos etal., 2008). However, in a structure thought 
to represent the closed CFTR state, the fourth intracel- 
lular loop is not in proximity of the residue correspond- 
ing to F508 (Hohl et al., 2012). A loss of stabilizing 
interactions present in the opening transition state, but 
not in the closed ground state, could explain why the 
major gating defect caused by disease mutation AF508 
is an increase in AG + ci— >oi, reflected by a dramatic 
(^30-fold; Miki et al., 2010) reduction in channel open- 
ing rate, and why P„ of this mutant is robustly stimulated 
by NPPB (Wang et al., 2005). 

NPPB slows channel closing by inhibiting the 
hydrolysis step 

The second effect of NPPB on WT CFTR gating was an 
approximately fourfold slowing of channel closing rate 
(Figs. 3, A and C, and 6, C and F) . Using maximum like- 
lihood fitting of burst duration distributions, we show 
that NPPB delays closure by slowing the Oi— >0? transi- 
tion (Fig. 5, C-E). Thus, NPPB bound at its "gating site" 
slows ATP hydrolysis at site 2 of the NBD dimer inter- 
face. This allosteric action implies a significant con- 
formational change at that NPPB binding site upon 
ATP hydrolysis at site 2, consistent with the proposal by 
Gunderson and Kopito (1995) of a TMD conforma- 
tional change associated with this step. Because the (likely 
very slow) rate of reversal of the Oi— >0? step cannot be 
measured (Csanady et al., 2010) , we cannot tell whether 
NPPB causes an isolated destabilization of the Oi— 
transition state (a pure "anticatalyst" effect), or whether 
the stability of the O2 ground state also changes re- 
lative to Oi. 
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Figure 1 2. Full kinetic model of NPPB effects on CFTR gating and permeation. (A; left) Two-tiered kinetic model of NPPB gating ef- 
fects (cube). States Ci, Oi, O2, and C2 (red) correspond to the four states depicted in the cartoons in Figs. 3-5 and represent states in 
which NPPB is not bound at its "gating site." The C2<->Ci transition, which reflects exchange of ADP for ATP at site 2, is modeled as a 
single step, with a of 50 pM for ATP. For simplicity, a single voltage-independent 1^ of 80 pM is used to characterize rapid binding/ 
unbinding of NPPB to the gating site in all four conformational states (vertical transitions; gray arrows). States Ci*, Oi*, Oj*, and C 2 * 
(blue) are conformational states analogous to Ci, Oi, O?, and C 2 , but with NPPB bound at the gating site. In 0 or 210 pM NPPB, the full 
model (cube) reduces to the four-state models to its left (red) and right (blue), respectively; states £ in the blue reduced model are 
compound states (C, = (Q; Ci*), 0 { = {Oi; Oi*), etc.), and printed rates are apparent rates of transition between them. (Right) Voltage 
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Prolongation of bursts by NPPB was already noted by 
Zhang et al. (2000) and interpreted as a "foot-in-the- 
door" mechanism: the presence of the blocker physically 
preventing gate closure. Several of our observations are 
inconsistent with such an assignment. First, unlike the 
pore block, the effect on burst duration shows little volt- 
age dependence (compare Figs. 3 C and 6 C). Second, a 
foot-in-the-door mechanism would be expected to slow 
the actual closing step, transition Og— >C%, for which we 
find no evidence (Fig. 5 E). Third, the presence of NPPB 
does not slow closure of the nonhydrolytic mutant 
K1250A (Fig. 4, A and B), confirming that the gate can 
close with NPPB bound in the pore, just as it can close in 
the presence of bound MOPS , as the latter does not af- 
fect the closing rate of either WT (Fig. 3, B and C) or 
K1250A CFTR (Fig. 4, A and B). Note that the similar 
effective valences of NPPB and MOPS" (Fig. 1, D, F, 
and H) suggest that their blocking sites are close to each 
other. The fact that the gate can close with either blocker 
bound is consistent with an extracellularly localized gate 
(Bai et al., 2011; Norimatsu et al., 2012), as predicted 
by the ABC transporter analogy, based on which CFTR's 
closed-pore conformation corresponds to inward-facing 
TMDs. Collectively, these results rule out a classical open- 
channel block mechanism as an explanation for the de- 
layed pore closure in NPPB. Moreover, the differential 
effect of NPPB on hydrolytic versus nonhydrolytic clos- 
ing rate (compare Figs. 3 C and 4 B) adds further sup- 
port for an underlying nonequilibrium gating cycle. 

The potentiator compound Vx-770, now approved for 
the treatment of CF patients carrying the G551D muta- 
tion, was proposed to slow closure of WT CFTR by facili- 
tating a reentry pathway from state Og to state Oi, which 
would allow ADP in site 2 to be replaced by a new 
ATP molecule without intervening pore closure (Jih and 
Hwang, 2013) . Such a mechanism would prolong steady- 
state burst durations by allowing occlusion and hydrolysis 
of more than one ATP molecule in site 2 within a single 
burst. In contrast, the macroscopic relaxation time con- 
stant after ATP removal measures mean burst duration in 
the absence of bath ATP, i.e., without the possibility of 
ATP rebinding. Because the macroscopic closing time 
constant in the presence of 210 pM NPPB (Fig. 6 C, 
blue time constant; T re i axation = 847 ± 60 ms; n= 14) is in- 
distinguishable from the steady-state mean burst duration 



(Fig. 6 F; T burst = 825 ± 145 ms; n = 15), we conclude that 
NPPB prolongs open bursts without significantiy chang- 
ing the near 1:1 coupling between ATP occlusion events 
and channel bursts of WT CFTR. 

Proposed NPPB mechanism of action can account for very 
potent stimulation of AF508 CFTR 

Our kinetic analysis accounts for the three- to fourfold 
stimulation by NPPB of P Q for WT CFTR, but does it 
account for the 10-15-fold stimulation of AF508 CFTR 
reported by Wang et al. (2005) — also reproduced here 
(Fig. 7) — which makes NPPB the most potent drug tar- 
geting this disease mutant available to date? Because P Q 
is a nonlinear function of mean burst and interburst du- 
rations (P (> = T b /(T b + Tjb)), such discrepancies between 
overall stimulation efficiencies and effects on micro- 
scopic rates are actually expected. Whenever the "con- 
trol" (without NPPB) value of P Q is low (T b « T ib ), the 
approximately fourfold increase in T b (Fig. 3 A; approxi- 
mately threefold at +60 mV; Fig. 6 F) and approximately 
fourfold decrease in T ib (Fig. 3 F; approximately three- 
fold at +60 mV; Fig. 6 F) should affect P (> independently, 
predicting an ^16-fold stimulation (approximately nine- 
fold at +60 mV) , in very good agreement with our obser- 
vations on AF508 CFTR (Fig. 7). This same argument 
also explains the observed dependence of stimulatory 
efficiency on phosphorylation level (Fig. 8 C; compare 
to Wang et al., 2005). Indeed, assuming a constant 
(entirely phosphorylation-independent) fourfold effect 
of 210 pM NPPB on both T b and T ib predicts ^15-, ~2-, 
and ^4-fold effects on P Q for unphosphorylated (before 
exposure to PKA, P„ of M).003), fully phosphorylated 
(in the presence of PKA, P Q of M).4) , and partially de- 
phosphorylated (after PKA removal, P„ of M).2) WT 
CFTR (see table in Fig. 8 D), in good agreement with 
our experiments (Fig. 8, A-C) . Note that at +60 mV, the 
approximately twofold increase in P Q of fully phosphor- 
ylated channels is largely counteracted by pore block, 
yielding little net effect (Fig. 8 B; compare to Wang 
etal., 2005). 

Simple unified kinetic model captures most key features 
of observed NPPB effects 

Our conclusions on the gating effects of NPPB can be 
summarized in a simple model that treats gating and 



and dose dependence of pore block by NPPB is modeled as an instantaneous effect on apparent unitary conductance (g) and is approxi- 
mated by the Boltzmann equation with parameters printed (go, control unitary conductance; T, temperature in Kelvin; R = 8.31 J • moP 1 • K _1 ; 
F = 96,500 C • moP 1 ). (B-L) Predictions of the model in A for the experimental protocols and analysis results obtained in this study. 
Macroscopic current time courses for the experimental protocols in B-E and I-L were calculated, whereas single-channel traces for F-H 
were simulated using standard Q-matrix techniques (Colquhoun and Sigworth, 1995). Rate Oi*— >0 2 * was tentatively set to zero. Gating 
effects of the K1250A mutation were modeled by setting rate Oi— >Oa to zero while increasing the Kj for ATP to 5 mM (Vergani et al., 
2003), those of the AF508 mutation were modeled by decreasing rate C[— >Oi 30-fold while increasing rate G^— >C] threefold (Miki et al., 
2010; Jih et al., 201 1). Analysis of macroscopic currents was done as described in Figs. 2-4, 6, 7, and 10. For F and H, five independent 
events lists, containing 100 open events each, were simulated for each condition; bar charts show mean ± SEM of obtained mean open 
and closed times. Open duration histograms in G were created from 3,000 simulated open events for both conditions and fitted by maxi- 
mum likelihood as described for Fig. 5 (C-E). 
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permeation effects separately (Fig. 12 A, left and right) . 
Gating (Fig. 12 A, left) is described by a simple two- 
tiered mechanism (Fig. 12 A, cube), in which the bot- 
tom tier represents gating states with NPPB not bound 
at the gating site and is identical to the four-state model 
depicted throughout the cartoons in Figs. 3-5. The 
C 2 — >Ci transition reflects the exchange of ADP for ATP 
at site 2, which in reality takes place in two sequential 
steps (ADP unbinding followed by ATP binding) . Be- 
cause under our recording conditions, in which the cy- 
tosolic face of the patch is continuously perfused with 
an ADP-free bath solution, unbinding of ADP is irre- 
versible (and likely rapid) , ADP-ATP exchange is mod- 
eled here as a single step, with a Kj for ATP of 50 uM 
(Csanady et al., 2000). Channel opening rate (rate 
Q— >0{) represents the inverse of the observed mean 
interburst duration (~1 s for WT CFTR; Figs. 5 B and 6 F) . 
Rates Oj— »0 2 and 0 2 — >C 2 reflect the fast and slow rate, 
respectively, obtained from fits to the distributions of 
open burst durations of WT CFTR (k Y and k 2 in Fig. 5 C) . 
Finally, rate C^— »Ci represents the slow rate of nonhy- 
drolytic closure, modeled by the closing rates of nonhy- 
drolytic mutants K1250A (Fig. 4 A) or El 37 IS (Fig. 1 K) , 
or of WT channels that have been locked open by ATP 
plus pyrophosphate (Fig. S2) ; the time constant of this 
slow process is ^30 s. 

Vertical transitions (Fig. 12 A, gray arrows) reflect 
binding/unbinding of NPPB to the gating site. Al- 
though our data do not provide direct estimates of these 
rates, the onset of the gating effects both upon the ad- 
dition and removal of NPPB appeared instantaneous 
(e.g., Figs. 3 D, 6 D, and 7, A and B) , without observable 
delay (which would have caused relaxation time courses 
to become sigmoidal) . This suggests that NPPB is at 
rapid equilibrium with the gating site, at least compared 
with the rate of our solution exchange, and hence rela- 
tive to the rates of CFTR gating. For simplicity, we as- 
sumed a single voltage-independent Kj of 80 uM in all 
four conformational states. 

States d*, Oj*, 0 2 *, and C 2 * of the top tier (Fig. 12 A, 
blue) are conformational states analogous to C 1; Oj, 0 2 , 
and C 2 , but with NPPB bound at the gating site. In the 
absence of NPPB, the full model (Fig. 12 A, cube) reduces 
to the four-state model to its left (red) . Moreover, because 
of the rapid equilibrium of NPPB with its gating site, in the 
presence of NPPB, each state of the bottom tier is at rapid 
equilibrium with the analogous state of the top tier, so that 
each pair can be treated kinetically as a single compound 
state. Thus, even in the presence of NPPB, the full scheme 
readily reduces to a four-state model with compound states 
C, = {d; Cj*}, 6, = {Oj; Oj*}, 6 2 = {0 2 ; 0 2 *|, and C 2 = (C 2 ; 
C 2 *(. Apparent transition rate X; -> Yj is obtained as 
//./.• • //"/(,*, where p- = p(x\%}, p* = p(x*|Xi), &y is rate 
X— >Yj, and Aj* is rate X* -» Y*\ fractional occupancies pi 
and p" are given by p x = K^/iK^ + [NPPB]) and p' = 
[NPPB] / (Aa + [NPPB] ) . In 210 uM NPPB, the full scheme 



(Fig. 12 A, cube) reduces to the four-state model to its 
right (blue): apparent rate C\ ->0! reflects the observed 
opening rate in 210 uM NPPB (Figs. 3 F and 6, B and F), 
rate 6[ -> 0 2 is the rate-limiting step for closure and re- 
flects closing rate in 210 pM NPPB (Figs. 3 C and 6, C 
and F) , 6 2 -> C 2 seemed unaffected by NPPB (see fits to 
the distributions of open burst durations in 20 pM NPPB; 
Fig. 5 D), whereas rate O, ->C, reflects K1250A closing 
rate in 210 pM NPPB (Fig. 4 B). The rates in the top tier 
of the full scheme were chosen to yield the above com- 
pound rates in the presence of 210 pM NPPB. 

Permeation effects of NPPB (Fig. 12 A, right) are mod- 
eled as simple voltage- and dose-dependent pore block 
and assumed instantaneous compared with gating effects. 
Thus, apparent unitary conductance (g) is approximated 
by the Boltzmann equation with parameters printed (g h 
control unitary conductance; T, temperature in Kelvin; 
R= 8.31 J • mof 1 • K" 1 ; F = 96,500 C • mol" 1 ; .M0) and z 
were chosen to best match fractional block measured at 
— 120 and +60 mV, at which the majority of our dataset 
was obtained) . 

The predictions of this simplified model were verified 
by calculating macroscopic current time courses cor- 
responding to all the experimental protocols tested in 
Figs. 2-4 and 6 and 7, as well as by simulating single- 
channel dwell-time sequences for all the conditions 
tested in Figs. 5 and 6 F. Macroscopic and single-chan- 
nel traces were analyzed exactly as the corresponding 
experimental traces, and results are summarized in 
Fig. 1 2 (B-L) . The predictions of die model were found to 
be in excellent agreement with our experimental obser- 
vations. Indeed, the model faithfully reproduced (a) 
robust stimulation of WT P„ as reflected by a discrep- 
ancy between fractional effects on steady-state macro- 
scopic and unitary currents (Fig. 12, B and C; compare 
to Figs. 2, A, B, and E, and 6, A and B) and current over- 
shoots upon rapid removal of NPPB (Fig. 12 E; com- 
pare to Figs. 3 D and 6 D) ; (b) slowing of WT (hydrolytic) 
macroscopic closing rate (Fig. 12 D; compare to Figs. 3, 
A and C, and 6 C); (c) acceleration of WT macroscopic 
opening rate (Fig. 12 E; compare to Figs. 3 F and 6 D); 
(d) shortening and prolongation, respectively, of steady- 
state single-channel mean closed (interburst) and open 
(burst) durations (Fig. 12, F and H; compare to Figs. 5 B 
and 6 F), the latter being caused by a slowed rate 
6, -^6 2 (Fig. 12 G; compare to Fig. 5, C and D); (e) ac- 
celerated closing rate of the nonhydrolytic K1250A mu- 
tant (Fig. 12 I; compare to Figs. 4, A and B, and 6 E), 
but (f) lack of effect on the P G of this channel (Fig. 12 J; 
compare to Fig. 4, C and D) , (g) observed apparent af- 
finities of the NPPB gating effects (Fig. 12 K; compare 
to Fig. 10) , as well as (h) extremely efficient stimulation 
of low P Q mutants such as AF508 (Fig. 12 L; compare 
to Fig. 7; also compare to Fig. 8). Thus, our simplified 
model captures most key features of the complex effects 
of NPPB on CFTR function. Indeed, the only feature 
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the model does not account for is the observed slight 
voltage dependence of the stimulatory efficiency: 210 uM 
NPPB affected both opening and closing rates by ap- 
proximately fourfold at —120 mV, but only by approxi- 
mately threefold at +60 mV. This small difference is 
little apparent when the control (unstimulated) P () is 
high. However, for low starting P„ values, NPPB-induced 
changes in the opening and closing rate affect P (> in a 
multiplicative manner, predicting an —16-fold stimula- 
tion at —120 mV, but only —9-fold at +60 mV, as found 
for the low P„ AF508 mutant (Fig. 7 D). This feature 
could be accounted for by assuming a slight voltage de- 
pendence of NPPB binding at the gating site; for in- 
stance, could be —80 uM at —120 mV but closer to 
-100 pM at+60mV. 

In conclusion, our analysis has identified two strategic 
intervention points in the CFTR gating cycle to robustly 
stimulate P„: stabilizing the transition state for opening, 
or slowing the ATP hydrolysis step. Both strategies 
uniquely rely on a nonequilibrium gating cycle. We 
have further identified a compound that is capable of 
exerting both effects. Identifying the NPPB gating site 
would be an exciting further step toward rational devel- 
opment of new nonblocking NPPB derivatives that stim- 
ulate CFTR gating with improved affinity, specificity, 
and/or potency. Such drugs — in combination with a 
suitable "corrector" compound to mend the folding de- 
fect and thermal instability of AF508 CFTR (Wang et al., 
2011; Rabeh et al., 2012) — could provide pharmacolog- 
ical treatment for >90% of CF patients carrying at least 
one AF508 allele. 
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